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. ^ , Abstract. We present data from high- resolution numerical simulations of the Navier- 

£^ I Stokes-a and the Leray-a models for two-dimensional turbulence. It was shown 

(— I ' previously (Lunasin et al., J. Turbulence, 8, (2007), 751-778), that for wavenumbers 

fS i' k such that ka ^ I, the energy spectrum of the smoothed velocity field for the two- 

dimensional Navier-Stokes-a (NS-a) model scales as k~'^ . This result is in agreement 
with the scaling deduced by dimensional analysis of the flux of the conserved enstrophy 
K^ I using its characteristic time scale. We therefore hypothesize that the spectral scaling 

of any a-model in the sub-a spatial scales must depend only on the characteristic time 
QQ , scale and dynamics of the dominant cascading quantity in that regime of scales. The 

O^l ' data presented here, from simulations of the two-dimensional Leray-a model, confirm 

our hypothesis. We show that for ka 3> 1, the energy spectrum for the two-dimensional 
Leray-a scales as k~^, as expected by the characteristic time scale for the flux of the 
conserved enstrophy of the Leray-a model. These results lead to our conclusion that 
the dominant directly cascading quantity of the model equations must determine the 
scaling of the energy spectrum. 
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1. Introduction 

In [20] we observed that the scahng exponent of the energy spectrum of the two- 
dimensional (2d) Navier-Stokes-a model (NS-a), for wavenumbers k such that ka ^ 1, is 
k'"^ . A posteriori, we saw that this scaling corresponding to that predicted by assuming 
that the dynamics for ka ^ 1 was governed by the characteristic time scale for flux 
of the conserved enstrophy. We were therefore led to speculate that (in general) the 
unknown scaling exponent for any a-model may be predicted by the dynamical time 
scales for the dominant conserved quantity for that model in the regime ka ^ 1. In this 
paper, we present new numerical simulations of the 2d Leray-a model which support 
this hypothesis. 

We measure the scaling of the energy spectra from simulations of two-dimensional 
flow, performed at a resolution of 4096^, in the limit as a — > oo, for two models: the 
NS-a model [10, 11, 15, 21] 

dtv — uAv — u xV X V = — Vp + /, (1) 

V ■u = V ■v = 0, (2) 
v = {I- a'^A)u, (3) 

and the Leray-a model [5] 

dtV - uAv + (n ■ V)v = -Vp + /, (4) 

V ■ M = V ■ t; = 0, (5) 
v = {I- a^A)u, (6) 

where v, u and p are the unsmoothed velocity, smoothed velocity and the pressure 
respectively for the Leray-a model and we use ~ to distinguish the variables in the NS-a 
model; v is the viscosity, and / is the body force. Notice that the two systems above 
reduce to Navier-Stokes equations (NSE) when a = 0. One can think of the parameter 
a as the length scale associated with the width of the filter which smooths v (or -li) to 
obtain u (or u). The filter is associated with the Green's function (Bessel potential) 
of the Helmholtz operator (J — a^A). We supplement both of the systems above with 
periodic boundary conditions in a basic box [0, L]^. 

The inviscid and unforced version of the three-dimensional (3d) NS-a was 
introduced in [15] based on the Hamilton variational principle subject to the 
incompressibility constraint div v = Q. By adding the viscous term —vAv and 
the forcing / in an ad hoc fashion, the authors in [1, 2, 3] and [10] obtain the 
NS-a system which they named, at the time, the viscous Camassa-Holm equations 
(VCHE), also known as the Lagrangian averaged Navier-Stokes-a model (LANS-a). In 
references [1, 2, 3] it was found that the analytical steady state solutions for the 3d 
NS-a model compared well with averaged experimental data from turbulent flows in 
channels and pipes for wide range of large Reynolds numbers. It was this fact which led 
the authors of [1, 2, 3] to suggest that the NS-a model be used as a closure model for 
the Reynolds averaged equations. Since then, it has been found that there is in fact a 
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whole family of 'a'- models which provide similar successful comparison with empirical 
data - among these are the Clark-a model [6, 7], the Leray-a model [5], the modified 
Leray-a model [16] and the simplified Bardina model [8, 19] (see also [23] for a family 
of similar models). 

The 3d NS-a model was tested numerically in [2] and [4], for moderate Reynolds 
number in a simulation of size 256^, with periodic boundary conditions. It was observed 
that the large scale features of a turbulent flow were indeed captured and there was 
a roll over of the energy spectrum from k~^^^ for ka <^ 1 to something steeper for 
ka ^ 1, although the scaling ranges were insufficient to enable extraction of the power 
law unambiguously. Other numerical tests of the NS-a model were performed in [13], 
[14], and [21], with similar results. 

In the limit as a — ^ oo, we call the two equations NS-oo [17] and Leray-oo, 
respectively, where the forcing term on both models are rescaled appropriately to avoid 
trivial dynamics. The equations for the NS-oo and Leray-oo are exactly the equations 
(1) and (4) together with the incompressibility condition except that equations (3) and 
(6) are replaced by the equations v = —L'^Au and v = —L'^Au, respectively. Under the 
assumption that the scaling of the spectrum as a — i> oo is identical to the scaling in the 
range ka ^ 1 for finite (small) a and sufficiently long scaling ranges, we obtain a high 
resolution numerical calculation of the sub-a scales. This assumption was verified in 
the finite a calculation of the NS-a model for two- dimensions in [20]. We stress again 
here, that one has to rescale the forcing appropriately in order to avoid trivial dynamics 
for large values of a, and decaying turbulence at the limit when a— ^oo [17]. 

Let Uk and V^ denote the typical smoothed and unsmoothed velocities of an eddy 
of size 1/k for the Leray-a model. Similarly, let Uk and Vk denote the typical smoothed 
and unsmoothed velocities of an eddy of size 1/k for the NS-a model. Such 'typical' 
velocities may be defined by the energy per unit area in the shell [k, 2k) as we will show 
in the next section. From our simulations of 2d NS-oo in [20], the energy spectrum of the 
smoothed velocity u scales as k'"^ . In this paper we will show from numerical simulations 
that the energy spectrum of the smoothed velocity u of the 2d Leray-oo model scales as 
k~^. These scalings can also be derived analytically (see [20] and section 2 below) under 
the assumption that an eddy of size 1/k, for ka ^ 1, has a typical time scale comparable 
to the inverse of the square root of the enstrophy contained in this eddy. That is, the 
dominant direct cascading quantity, which is the enstrophy in the 2d NS-a and Leray- 
a, dictates these typical time scales. Specifically, under this assumption, the governing 
time scales for an eddy of size 1/k in each model are given by {kVk)~^ (for NS-a) and 
{k\/UkVk )~^ (for Leray-a). We assert that the difference in the dominant forward 
cascading conserved quantities in these two models is what leads to the different power 
laws. Our numerical results in 2d for two different a-models, with different forward 
cascading conserved quantities, support this assertion. 

Based on our studies in 2d, we extrapolate our conclusions to the 3d case as follows. 
For the 3d NS-a model and 3d Leray-a model, the governing time scales for an eddy of 
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size 1/k, for ka ^ 1, must be given by {kyUkVk )~^ and {kVk)~^, respectively. This is 
because the energy conserved (in the absence of forcing and viscosity) in the 3d NS-a 
is given by J^u -v dx while in the 3d Leray-a it is given by J^v -v dx. Accordingly, we 
assert that for /ca ^ 1, the energy spectra of the smoothed velocity fields in the 3d case 
will scale as k~^^^^ (steeper than k~^ as originally suggested in [11]) for the 3d NS-a, and 
as /c~^^/^ (steeper than /c~^^/^ proposed in [5]) for the 3d Leray-a model. This assertion 
is yet to be confirmed computationally in future work. Our prediction of /c~^^/^ power 
law for the smoothed energy spectrum of the 3d Leray-a model corresponds to one of 
the three candidate power laws derived in [5] . The idea that the average velocity of an 
eddy of size of the order 1/k can be evaluated in three different ways, which will then 
lead to three different power laws, was in fact first introduced in [5]. 

Throughout the paper we denote by r^ the characteristic time scale of an eddy of 
size 1/k, Q = [0,L]^. We denote the rough and smoothed vorticities by q = V x v 
(or q = V X v) and u = 'V x u (or u) = V x -u), respectively. The paper is organized as 
follows. In section 2 we derive the power laws for the 2d Leray-a model and then give 
a comparison to the corresponding power laws of the 2d NS-a equations. In section 3, 
we give a brief review of the numerical results in [20] and then present our numerical 
results for the 2d Leray-a. In the last section, we give a summary of our main results 
and give a brief description of how this study can help us predict the unknown power 
laws for the 3d NS-a and 3d Leray-a equations. As mentioned above, our predictions, 
based on this study, on the power laws for the two models just mentioned, are different 
from those suggested in [5, 11]. In those works, the k~^ and k~^^'^ power laws for the 
3d NS-a and 3d Leray-a model, respectively, were proposed under the assumption that 
the time scale which governs the small scales is determined by the smoothed velocity 
field alone (even though there were two other candidate power laws derived in [5]). 

We dedicate this work to our friend and colleague Darryl D. Holm on the occasion 
of his 60th birthday in acknowledgement of his continuing support and inspiration in 
stimulating scientific interactions and discussions over the past years and many to come. 

2. Navier-Stokes-a vs. Leray-a model in two dimensions 

In this section we give a comparison between the two a-models. For completeness we 
briefly present the analytical arguments for the different power laws of the energy spectra 
which arise in the 2d Leray-a equations. For complete details we direct the reader to 
look at the derivation of power laws of the 2d NS-a in [20] (see also [5, 6, 8, 9, 11, 12, 16] 
for the analytical calculation of the power laws of energy spectra for the other a-models). 
To compute the scaling of the smoothed energy spectrum E'^{k) of the 2d Leray-a in 
the wavenumber regime ka '^ 1 in the forward enstrophy inertial subrange, we start 
by splitting the flow into the three wavenumber ranges [1, /c), [A;, 2A;), [2/c, oo). For a 
wavenumber /c, we define the component Uk of a velocity field u by 

Uk := uu{x) = ^ w(Oe*^«-^ (7) 
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and the component Uk'.k" with a range of wavenumbers [k', k") by 

Uk',k" ■=Uk',k"ix) = ^ Uk. (8) 

k'<k<k" 

E'^{k) is then the energy spectrum associated to 

l/ll l|2 

^fc - 2 \\\'^k,2k\\L2{Q) 

which is the average (with respect to an infinite time average measure [12]) smoothed 
energy per unit mass of eddies of hnear size / G (^, ^]- 

We assume kf < k, where kf is the forcing wavenumber, since we are interested 
on the effects of the Leray-a model in the enstrophy cascade regime. We decompose 
the u, V and the V x u and V x t> corresponding to the three wavenumber ranges. 
We then write the enstrophy balance equation for the Leray-a model for an eddy of 
size ~ k~^. Taking an ensemble average (with respect to infinite time average measure) 
of the enstrophy balance equation we get 

u [{\\Auk,2k\\hin) + "1 VAnfc,2fc||i2(^))) = (Z,) - (Z^k) , (9) 

where Zk may be interpreted as the net amount of enstrophy per unit time that is 
transferred into wavenumbers larger than or equal to k. Similarly, Z2k represents the 
net amount of enstrophy per unit time that is transferred into wavenumbers larger than 
or equal to 2k. Thus, Zk — Z2k represents the net amount of enstrophy per unit time that 
is transferred into wavenumbers in the interval [k.,2k). We then rewrite the averaged 
enstrophy transfer equation (9) as 

ueE^{k) ~ u / k^E^{k)dk ~ {Zk) - {Z^k) , 



where E^ik) is the energy spectrum associated to 

^k = 2 \ll"fc.2fc|li2(f^) + a^||VMfc,2fc||i2(n)) , (10) 

which is the average energy per unit mass of eddies of linear size / G ij^,^]- 

Thus, as long as uk^Ea{k) <^ {Zk) (that is, {Z2k) ~ {Zk)i there is no leakage of 
enstrophy due to dissipation), the wavenumber k belongs to the inertial range. In the 
forward cascade inertial subrange, we follow Kraichnan [18] (see also [9]) and postulates 
that the eddies of size larger than 1/k transfer their energy to eddies of size smaller 
than l/{2k) in the time r^ it takes to travel their length ~ 1/k. That is, 

kUk 
where Uk is the average velocity of eddies of size ~ 1/k. Since there are two different 
velocities in the Leray-a model, there are three physically relevant possibilities for this 
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U', 



-^ / \vkMx)\^dx\ ~ ( / (1 + a^e)E^{k)dk] ~ (A;(l + a^e)E^{k)y 

/ I f \ 1/2 / /.2fe \ 1/2 

Ul = (j^ / "fc.2fc(a;) ■Vk,2kix)dx) ~ ( / Ea{k)dk\ ~ {kEa{k)f^'^, 



.2 _ / 1 r ., ..^,2. \^^^ . . / r i^«(^) .A ' . . ( kE^ik) 



Ut = (^ \UkM^)\'dx ~ ,. . \ o^ ^fe 



These define the aforementioned 'typical' velocities, in particular Vk = U^ and Uk = U^. 
Corresponding definitions may be made for Vk and Uk using the variables for the NS-a 
model, see [20]. Thus, 

We may therefore write the typical time scale of an eddy of size k~^ in (11) as 

That is, in keeping with the historical approach to the problem [9, 10], there are in 
principal three different time scales and it is left to empirical evidence to infer the 
correct time scale for a particular a-model. 

The enstrophy dissipation rate rja which is a constant equal to the flux of enstrophy 
from wavenumber k to 2k is given by 



4i ^^^-w^-^ikif^- ("' 

r/^/'(l + a2fc2)(n-l)/3 



and hence 

EM ^3 

Thus, the energy spectrum of the smoothed velocity u is given by 



F (k\ { ^c/^ ' ^^^^ to < 1 , 

E^{k)= "^, , ^ <^ 2/3 (15) 

^' l + a^F ] ^^^-(i7-2n)/3 when to »1. ^ ' 

Therefore, depending on the average velocity of an eddy of size k~^ for the Leray-a 
model, we obtain three possible scalings of the energy spectrum, /^-{i7-2n)/3^ ^^ _ g^ ^^ 2) 
all of which decay steeper than the Kraichnan k~^ power law, in the subrange ka ^ 1. 
The goal, which we stress here again, is to infer the correct time scale by measuring the 
scaling exponent of energy spectra computed from high-resolution numerical simulation 
data. 

We summarize some points of comparison between the two models in Table 1. In the 
absence of viscosity u and the forcing /, the two conserved quantities, namely the energy 
and enstrophy, for the two models, in two dimensions, are specified in the first block- 
row of Table 1. Since the energy in 2d flow goes upscale ([9, 12, 18, 20]), we are more 
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interested in the enstrophy which has its dominate cascade downscale ([9, 12, 18, 20]). 
For the NS-a model, the enstrophy conserved is Qa = Jq \(j\'^dx, while for the Leray-a, 
the conserved enstrophy is given hj Q^ = J^q ■ u dx. It is for this reason that the 
characteristic time scale for eddies of size smaller than the length scale a, for the two 
models, differ. The second block-row gives the three possible characteristic time scales, 
and the corresponding scaling predictions, three for each of the models. The notation 
in Table 1 defaults to that for NS-a; replace Uk and Vk by Uk and Vk, respectively, in 
the formula for Tk for Leray-a. 

These three possibilities for t^ and the corresponding power laws are given in 
Table 1. In [20], it was shown that the energy spectrum of the 2d NS-oo equations 
attains a power law of k~'^ as the resolution is increased. The convergence of the 
spectral scaling is presented in the third block-row in Table 1. A similar study for 
the 2d Leray-oo shows a convergence to a power law of k~^^ as we shall show in the next 
section. 



Table 1. Comparison between the NS-a and Leray-a in two dimensions (in the second 
block-row, Uk and Vk should be replaced by Uk and Vk respectively to obtain the Tk 
for Leray-a). 

NS-a Leray-a 

Ideal invariants: Energy e^, = J„ u • w dx Energy e^ = J^u ■ v dx 

Enstrophy fta = /jj \q\'^ dx Enstrophy fi^ = J^q ■ oj dx 
Expected scaling 
in the range ka ^ 1 

fc-5 
fc-13/3 



5.5 
5.2 
5.0 



ifTfc = (fcl4)-l 


fc-7 


ifTfc- (kJUkVk 

V / 


fc-19/3 


ifTfe = (fct7fe)-l 


fc-17/3 


Convergence of k''^ 




as resolution is increased 




1024^ 7 = 


7.4 


2048^ 7 = 


7.1 


40962 7 ^ 


7.0 



3. Numerical results 

3.1. Details of the numerical simulation 

The Leray-a equations were solved numerically in a periodic domain of length L = 1 on 
each side. The wavenumbers k are thus integer multiples of 27r. A pseudospectral code 
was used with fourth-order Runge-Kutta time-integration. Simulations were carried out 
with resolutions ranging from 1024^ up to 4096^ on the Advanced Scientific Computing 
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QSC machine at the Los Alamos National Laboratory. To maximize the enstrophy 
inertial subrange, the forcing is applied in the wavenumber shells 2 < k < 4. We 
also add a hypoviscous term n{—A)~^v which provides a sink in the low wavenumbers. 
To discern a clear power-law of the Leray-a model spectrum, we consider data from 
simulation of the Leray-cxo equations 

dtv - {u ■ V)v = -Vp + uAv + f (16) 

V-n = V-w =0 (17) 

V = -L'^Au, (18) 

Similar to the case of the 2d NS-oo equations [17, 20], this allows us to see the scaling 
of the Leray-a model energy spectrum without contamination by finite-a effects. 

3.2. Results for the Leray-oo equations 

In Figures 1 and 2 we use the notation E{k) for E^{k). In Figure 1 we present the 
main numerical results from [20] for the 2d NS-oo equations, showing the A;~^ scaling 
of the energy spectrum. In Figure 2 we show that the 2d Leray-oo energy spectrum 
attains a k~^ power law in a 4096^ resolution simulation. From Table 1, the scaling k~^ 
stems, based on the analytical arguments in section 2, from a characteristic time scale 
given by r^ = {ky/VkUk) , which is comparable to the inverse of the square root of the 
enstrophy of an eddy of the size 1/k. To see this, recall from section 2 that the typical 
smoothed and unsmoothed velocity of eddy of size ~ 1/k are given by 

Uk = J^{\\Uk,2ki^)\\L2{Q)) , (19) 

Vk=^{\\Vk,2k{^)\\L2^n)), (20) 

then, we can define 

^k = j^ iW"^ ^ Uk,2k\\LHn)) , (21) 

Qfc = -^(||V X Wfc,2fclU2{Q)>, (22) 

as the smoothed and unsmoothed enstrophy per unit area in the shell [k,2k). Now, 
observe that 

rk= [k^/vdhY' -- {vokmy'. 

Therefore, the numerical results in Figure 2 clearly support our claim that the 
characteristic time scale t^ = (ky/VkUk) determined by the dominant cascading 
quantity, namely the enstrophy J^q ■ u dx for the 2d Leray-oo model, governs the 
dynamics of eddies in the subrange ka ^ 1. This conclusion is consistent with our 
original prediction in [20] that, in general, the dominant cascading conserved quantity 
in the a-model dictates the time scale associated with eddy of size 1/k and hence the 
power law of the energy spectrum E^{k) in the subrange ka^ 1. 
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Figure 1. Energy spectra E{k) = E"'{k) from a 4096^ simulation of the NS-a 
equations, (a): the blue line is the spectrum for NS-oo, the black line for the 
NSE (a = 0). (b) top to bottom: the energy spectrum for NS-oo compensated by 
/c^, k^^/^^ and fc-"^"^/^, respectively. The compensated spectrum in the top panel is flat in 
the range 6 < fc < 40, indicating the nominal range over which the fc~^ scaling holds. 
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Figure 2. Energy spectrum E{k) = E'"'{k) from a 4096^ simulation of the Leray-a 
equations, (a): energy spectrum for the Leray-oo. (b) top to bottom: energy spectrum 
of Leray-a compensated by fc-'^^/^, fc^, and fc^"^/^, respectively. The compensated 
spectrum in the middle panel is flat in the range 7 < A: < 50, indicating the nominal 
range over which the k^^ scaling holds. 



4. Conclusion 



The main goal of tliis study is to verify our claim in [20] about the choice of particular 
characteristic time scales of eddy of size 1/k, for /ca ^ 1, for particular a- model 
equations. Our results in [20] led us to conclude that the choice depends on the form of 
the cascading conserved enstrophy, which is the dominant forward cascading quantity. 
To verify this conclusion, we perform a high resolution simulation of the 2d Leray-a 
equations in the limit as a— ^oo similar to our study of the 2d NS-oo equations. We 
summarize the three steps to this study which verifies these claims: 
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(i) Identify the conserved quantities (in the absence of viscosity and forcing) for 2d 
Leray-a, and the dominant one in the forward cascade regime. 

(ii) Calculate the power laws using semi-rigorous arguments as in [6, 8, 9, 11, 12, 16, 20]. 
This will give us the three possible power laws of the energy spectrum in the 
wavenumber regime /ca 3> 1. 

(iii) Perform a high-resolution simulation to identify which one of the power laws of the 
energy spectrum calculated in step (ii) actually arise for the 2d Leray-a model. 

As we have speculated in [20], the scaling exponent in the wavenumber regime 
ka ^ 1 will be governed by the time scale of the dominant cascading conserved 
quantity in that regime. If we extend the same argument to predict the scaling 
for the 3d NS-a and the 3d Leray-a model equations then we obtain the following 
predictions. Since i?™ = | L ^aU-v dx and E'"'" = | L ^aV-v dx are the conserved 
energy which are the dominant cascading quantities for the 3d NS-a and 3d Leray- 
a, respectively, then we predict the scaling of E^{k) ~ k~^^'^ for the 3d NS-a (that 
is, steeper than the k~^ proposed in [11]) and the scaling of £""(A;) ~ k~^'^^^ for the 
3d Leray-a (that is, steeper than the k~^^^^ suggested in [5] in the wavenumber 
regime fca ^ 1. Our prediction of k~^'^^^ power law corresponds to one of the 
power laws initially derived in [5] as candidate power law for the smoothed energy 
spectrum of the 3d Leray-a model. The verification of these possibilities in the 3d 
case will be explored in future work. 
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